A Constraint on Electromagnetic Acceleration of Highest Energy Cosmic Rays 
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The energetics of electromagnetic acceleration of ultra-high-energy cosmic rays (UHECRs) is 
constrained both by confinement of a particle within an acceleration site and by radiative energy 
losses of the particle in the confining magnetic fields. We demonstrate that the detection of ~ 
3 x f 20 eV events is inconsistent with the hypothesis that compact cosmic accelerators with high 
magnetic fields can be the sources of UHECRs. This rules out the most popular candidates, namely 
spinning neutron stars, active galactic nuclei (AGNs), and 7-ray burst blast waves. Galaxy clusters 
and, perhaps, AGN radio lobes remain the only possible (although not very strong) candidates for 
UHECR acceleration sites. Our analysis places no limit on linear accelerators. With the data from 
the future Auger experiment one should be able to answer whether a conventional theory works or 
some new physics is required to explain the origin of UHECRs. 

PACS numbers: 41.60.-m, 96.40.-z 



I. INTRODUCTION 

The detection of ultra-high-energy cosmic rays (UHE- 
CRs) with energies above 10 20 eV has posed a challenge 
to the understanding of their origin and nature. At 
present, 17 such events were reported by the AGASA 
group and 2 events were observed by HiRes an energy of 
the one of them is estimated to be ~ 3 x 10 20 eV, close to 
the energy of the largest (~ 3.2 x 10 20 eV) event observed 
with the Fly 's Eye detector Nearly isotropic distribution 
on the sky and the absence of large-scale clustering sug- 
gests the cosmological origin of UHECRs. This is in ap- 
parent conflict with the observed energy spectrum which 
atsepi: 
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traveled the distance smaller than r- 
for protons with E > 3 x 10 20 eV). 

Cosmic rays are accelerated in astrophysical sources ei- 
ther by repeated scattering off macroscopic flows, such as 
shocks, winds and outflows, turbulent flows, or directly 
by an induced electric field around a magnetized rotat- 
ing object. For the acceleration to operate, a particle 
must remain confined within the acceleration region: the 
gyro-radius of the particle should not exceed the size of 
the system, R. This sets the maximum energy of the 
accelerated particle, 



lacks the Greisen-Zatsepin-Kuzmin (GZK) cutoff at en- 
ergy - 5 x 10 19 eV [HE! indicating that UHECRs have 

60 Mpc (~ 20 Mpc 



ZeBR ~ 9.3 x 10 23 Z B R kpc eV, (1) 



where Ze is the charge of a particle, B is the characteris- 
tic magnetic field strength in the acceleration region (in 
gauss), and i?kpc is the size in kiloparsecs. If the whole 
medium is moving relativistically with the Lorentz factor 
r towards an observer, E acc is boosted to r_B acc , with B 
and R being measured in the co-moving frame. Speci- 
fying the particle energy, say _E acc = 3 x 10 20 eV (the 
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largest observed), and assuming Z ~ 1, one obtains the 
magnetic field-size relation 0,0], which is plotted in Fig- 
ure^with the long-dashed line. Any astrophysical object 
to the right from this line can accelerate protons to en- 
ergies > 3 x 10 20 eV. 



Various energy losses, such as collisional and in- 
verse Compton scattering in a radiation field, limit the 
UHECR energy. But even in the absence of any such pro- 
cesses, an energetic particle will still loose energy radia- 
tively while moving through magnetic and electric fields 
which confine and accelerate this particle. The radiative 
(synchrotron) losses have been considered by various au- 
thors 0, 0, El, who usually estimated the synchrotron 
cooling time, but didn't take into account the finite size 
of the source self-consistentlv. fuij 

In this paper we re-analyze the energetics of electro- 
magnetic acceleration to elucidate its inherent and in- 
evitable limitations. The Hillas criterion, equation QJ, 
is shown to be often inaccurate and even misleading. We 
consider the most idealized models that carry only the 
most robust properties of cosmic accelerators, any de- 
tails of the electromagnetic acceleration process are not 
important for us. Thanks to the simplicity of the mod- 
els, analytical solutions for the particle energy, which ac- 
count for a source size self-consistently, were obtained. 
These results represent the most relaxed (and, hence, un- 
avoidable) constraints on the maximum energy of an elec- 
tromagnetically accelerated particle. This implies that 
the most favorable sources of UHECRs, such as neutron 
stars (NSs), active galactic nuclei (AGNs) and gamma- 
ray bursts (GRBs), cannot accelerate protons to the en- 
ergy ~ 3 x 10 20 eV and, hence, must be ruled out. How- 
ever, our critetion does not apply to linear accelerators 
(e.g., axial jets) 0. The Auger experiment may be ca- 
pable of ruling out the remaining candidates, radio lobes 
of AGNs and galaxy clusters, and if it does, the whole 
conventional astrophysical picture of UHECRs as accel- 
erated particles must be revisited. 
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II. INEFFICIENT ACCELERATION 

Let us consider diffusive acceleration first. This type of 
acceleration operates in shocks of 7-ray bursters, galaxy 
clusters, jets from AGNs interacting with the intergalac- 
tic medium and producing radio lobes, and, perhaps, in 
AGN cores near the base of a jet. Let us consider a 
shock propagating through a magnetized medium (inter- 
cluster medium [ICM], for instance), as shown in Figure 
|5Ji. Magnetic field may be inhomogeneous on a scale 
~ R. An accelerated particle gains energy by repeated 
scattering off a shock or a flow. After every scattering, 
the particle travels a great distance along the Larmor 
orbit until it returns and gets another kick. As long as 
the particle moves freely in the magnetic field it radiates 
and slows down. We will see that the maximum terminal 
energy of a particle (i.e., when the particle escapes the 
system), E, is determined by these radiative losses, but 
is insensitive to how large the energy, Eq, of the particle 
at the shock front is. Therefore, we refer to this regime 
as "inefficient acceleration". 

Let us consider a particle with some initial energy Eq 
propagating through a region of a size R filled with a 
magnetic field B. The energy of the particle gradually 
decreases according to the equation (see @): 



dE 

dx 



-FRad 



Ze 



B 2 {x)E\ (2) 



where x is the distance along the particle trajectory, -FR ac j 
denotes the radiation friction force, and Am p is the par- 
ticle mass. In @ only the transverse field component 
enters; we assimed that B± ~ B. The solution of this 
equation is 



E- 1 = V 



e: 



(3) 



lob = ZeB/Arripc), that is for the field strengths not ex- 
ceeding (Amp) 2 c 4 / (Z e) 3 in the rest frame of a particle. 
This yields the condition for the particle energy 

£« ( 1-1.9 x W 31 (A/ ZfB- 1 eV, (5) 

Ze J B 



which is satisfied for E < 10 21 eV for practically all 
sources. 



III. EFFICIENT ACCELERATION 

Let us now consider "cosmic inductors" where particles 
are accelerated by electric fields induced by rapid rotation 
of a magnetized object. Acceleration of this type should 
occur in neutron star magnetospheres and around accret- 
ing supermassive black holes in the centers of AGNs. We 
naturally assume that the magnetic field has a dipolar 
(or a multipolar) structure, hence field lines are bent on 
a scale ~ R, as shown in Figure |2t>- Because of rapid ro- 
tation, there is an induced electric field E in( i ~ |v x B|/c 
which accelerates a particle. The maximum value of -Emd 
is achieved near the light cylinder is close where it 
is equal to, at most, £"i n d ~ B. The particle in such a 
system gains energy rapidly, within one passage through 
the system. Hence one must retain the electromagnetic 
accelerating force Fem = Ze E- ln< i ~ Ze B in the energy 
equation |3J). Then it reads 



dx 3 \Am p c 2 



B 2 E 2 . 



(0) 



where Eq is the initial energy of the particle and 



For a small initial energy of an accelerated particle (Eq <C 
E acc , E cr ), the solution of this equation takes a simple 
and elegant form: 



E r 




B 2 (x) dx 



2.9 x 10 16 \' ' eV. 



B 2 R 



kpc 



For simplicity, it is assumed here that B(x) ~ constant 
within the system. It must be clear now that no mat- 
ter how energetic the particle is (Eq — * 00), after trav- 
eling through a region with a magnetic field its energy 
will not exceed the critical energy _E cr . Specifying the 
energy E CI = 3 x 10 20 eV and assuming A ~ 1, one ob- 
tains another B vs. R constraint, shown in Figure^with 
the short-dashed line. All astrophysical sources located 
above this line have E CI < 3 x I0 20 eV and hence cannot 
accelerate UHECRs. 

It should be mentioned that the classical expression 
@ for the radiation friction force is valid if the wave- 
length of the emitted radiation is larger than the "classi- 
cal radius" of a charge, c/lub S> (Z e) 2 / (Am p c 2 ) (where 



E — y E acc E cr tanhi/ E acc /E c 



(7) 



where E is the terminal energy of the particle. The so- 
(4) lution has two obvious asymptotics. If E acc ^ E cl - one 
recovers the Hillas constrain |4j E ~ E acc (equation P), 
whereas in the opposite limit one has 



E„ 



\j E acc E cr 

1.3 x IQ 20 A 2 Z~ 3/2 B- 1/2 eV. 



(8) 



This equation, in fact, follows from the balance between 
the acceleration and radiative losses, Fem = -Fkad, m 
equation JJ]J. Hence, we refer to this regime as "efficient 
acceleration" . Acceleration above this energy is impossi- 
ble because at larger energies radiative friction begins to 
dominate over electromagnetic acceleration. This con- 
straint for £ mM = 3 x 10 20 eV is plotted in Figure Q] 
with the dot-dashed line. No objects above this line can 
accelerate cosmic rays to this energy. 
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IV. DISCUSSION 

The conventional astronomical picture for the origin 
of these cosmic rays is the acceleration of charged par- 
ticles, e.g., protons or heavier atomic nuclei, in extra- 
galactic objects H- There are only few types of such 
objects. The most favorable are: spinning neutron stars 
(NSs) and magnetars central regions of active galactic 
nuclei (AGNs) AGN radio lobes and 7-ray burst (GRB) 
shocks. It is unlikely that shocks in galaxy clusters are 
the UHECR sources because they are too far, beyond 
the GZK distance, and they are not able to accelerate 
protons to energies above few times 10 19 eV [Toj| . 

In the previous sections we demonstrated that the en- 
ergy of a particle confined by magnetic fields within an 
acceleration site is determined by its radiative losses. We 
now discuss the constraints on B and R given by equa- 
tions CJ, Q, and (JHl and plotted in Figure ^ 

First, the constraint (JHJ shown with the dot-dashed 
line is the most stringent. It tells that protons cannot 
be accelerated electromagnetically at the sources which 
lie above this line. Equation JSJ also holds for heavier 
(e.g., iron) nuclei. Hence, compact stars, AGN cores, 
and GRB shocks (except during the late afterglow phase) 
are readily ruled out from the list of possible sources of 
ultra-high-energy protons and nuclei. 

Second, the remaining candidates, i.e., the radio lobes 
and galaxy clusters, may accelerate particles only via the 
diffusive mechanism, hence equation Q is appropriate. 
This equation together with specifies the allowed B— 
R parameter region. In the figure the dotted line corre- 
sponds to iron nuclei with the energy 3 x 10 20 eV and 
the solid lines correspond to protons with three energies: 
3xl0 20 , 10 22 , and3xl0 23 eV. One can see that the AGN 
radio lobes are at most marginally consistent with being 
the sources of the highest energy cosmic rays. Moreover, 
only a handful of such sources are relatively nearby (e.g., 
M87, Cen A, NCG 315) but their angular distribution is 
completely uncorrelated with the nearly isotropic distri- 
bution arrival directions of UHECRs. Shocks in galaxy 
clusters can be such sources, according to Figure ^ It 
has been argued from a more detailed analysis of shock 
acceleration, however, that they are not able to accel- 
erate protons above ~ 6 x 10 19 eV (13]. Overall, large 
objects are more preferable candidates for the highest en- 
ergy cosmic ray sources, both due to the larger terminal 



energy of an accelerated particle and the larger energy 
reservoir available, see Fig. 5 in the paper by Kronberg 

El 

Equation JTJ together with J3J or (jHJ puts a lower 
bound [16j on the size of the accelerating source: 

R > 6.0 x 1CT 5 £f Z 2 A- 4 kpc, (9) 

where E20 = E/10 20 eV. This is a quite remarkable re- 
sult since it sets the absolute limit on R, independent of 
the field strength. Now, two special cases follow. First, 
the size R exceeds the GZK distance ~ 20 Mpc when the 
energy of the proton is larger than -Egzk ~ 7 x 10 22 eV. 
That is, such an energetic proton will loose its energy 
through the interaction with the 2.7 K background radi- 
ation right at the acceleration site. Hence, it is unlikely 
that protons may be accelerated to the energies above 
Egzk- Second, above the energy i? H or ~ 4 x 10 23 eV, 
the size of the accelerator exceeds the size of the Uni- 
verse. Thus, i?Hor is the ultimate upper bound on the 
energy of electromagnetically accelerated protons. Note 
that this energy limit is below the energy of primary pro- 
tons (~ 10 24 eV) in the most popular (among other mod- 
els 

Z-burst model. 
To conclude, we arrived at an interesting result. Prac- 
tically all known astronomical sources are not able to pro- 
duce cosmic rays with energies near few times 10 20 eV. 
There is not too much room left for the conventional elec- 
tromagnetic (in a broad sense) acceleration. Pushing ob- 
servations up in energy by about an order of magnitude 
should clarify this situation. During ten years of opera- 
tion of the Auger cosmic ray observatoryabout 300 events 
above 10 20 eV are expected to be recorded. Extrapolat- 
ing the present AG AS A spectrum one can then expect a 
few events near or above 10 21 eV to be recorded, which 
may be enough to exclude radio lobes from UHECR 
sources. Should this happen, a considerable revision of 
the current astrophysical picture will be inevitable. 
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FIG. 1: The B vs. R diagram for UHECR sources. The 
long-dashed line is the original Hillas relation, eqn. for 
a proton of energy 3 x 10 20 eV. The short-dashed and dot- 
dashed lines represent the radiative cooling constraints for 
the diffusive and inductive acceleration, given by eqns. Ql 
and JSJ, for the same proton energy. Note that above the 
dot-dashed line the force of radiative friction dominates over 
any electromagnetic forces. The dotted and solid lines rep- 
resent the boundaries of the allowed parameter regions for 
3 x 10 20 eV iron nuclei and for protons of energies 3 x 10 20 , 



icr 



and 3 x 10 eV, respectively. Only those astronomi- 



cal objects which fall inside the "wedges" are, in principle, 
capable of accelerating the particles to such energies. The 
gray vertical lines mark two characteristic scales: the GZK 
attenuation distance (~ 20 Mpc) and the Hubble horizon size 
(~ 4 Gpc). For GRBs we took into account that the Lorentz 
boost changes with radius. 



FIG. 2: Cartoons representing a typical system with ineffi- 
cient (diffusive) acceleration (a) and with efficient (inductive) 
acceleration (b). 



